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Abstract 

Research on regulatory mechanisms in biological populations often focuses on environmental covariates. An integrated 
approach that combines environmental indices with organismal-level information can provide additional insight on 
regulatory mechanisms. Survival of spring/summer Snake River Chinook salmon {Oncorhynchus tshawytscha) is consistently 
low whereas some adjacent populations with similar life histories experience greater survival. It is not known if populations 
with differential survival respond similarly during early marine residence, a critical period in the life history. Ocean 
collections, genetic stock identification, and otolith analyses were combined to evaluate the growth-mortality and match- 
mismatch hypotheses during early marine residence of spring/summer Snake River Chinook salmon. Interannual variation in 
juvenile attributes, including size at marine entry and marine growth rate, was compared with estimates of survival and 
physical and biological metrics. Multiple linear regression and multi-model inference were used to evaluate the relative 
importance of biological and physical metrics in explaining interannual variation in survival. There was relatively weak 
support for the match-mismatch hypothesis and stronger evidence for the growth-mortality hypothesis. IVlarine growth and 
size at capture were strongly, positively related to survival, a finding similar to spring Chinook salmon from the IVlid-Upper 
Columbia River. In hindcast models, basin-scale indices (Pacific Decadal Oscillation (PDO) and the North Pacific Gyre 
Oscillation (NPGO)) and biological indices (juvenile salmon catch-per-unit-effort (CPUE) and a copepod community index 
(CCI)) accounted for substantial and similar portions of variation in survival for juvenile emigration years 1998-2008 (R^> 
0.70). However, in forecast models for emigration years 2009-2011, there was an increasing discrepancy between 
predictions based on the PDO (50-448% of observed value) compared with those based on the NPGO (68-212%) or 
biological indices (CPUE and CCI: 83-172%). Overall, the PDO index was remarkably informative in earlier years but other 
basin-scale and biological indices provided more accurate indications of survival in recent years. 
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Introduction 

One focus of population ecology is the identification of 
environmental indices that are related to variation in population 
size or productivity [1,2,3]. Such relationships are often based on 
hypothesized mechanisms, such as a "stable ocean" [4] or 
"optimal upwelling window" [5], but the relationships fail to hold 
up over time [6,7]. A parallel approach has been to identify a 
"critical period" in a species' life history, after which variation in 
the rate of mortality stabilizes [8,9]. If the critical period is 
successfully identified, then the abundance or condition of a cohort 
during, or shortly after, this critical period should provide a robust 
indication of relative survival. This approach is not necessarily 
based on a mechanistic understanding of mortality but can focus 
research efforts by identifying the critical period(s) in a species' life 



history. Moreover, a combined approach can identify relevant 
local or regional environmental factors and also provide insight on 
the timing and mechanisms of major mortality events. 

The Columbia River basin is the largest watershed on the west 
coast of the United States [10] and supports numerous populations 
of anadromous Chinook salmon, including five that are currendy 
listed under the U. S. Endangered Species Act (ESA) [11]. 
Extensive modifications have been made to the river's hydropower 
system to minimize mortality of juvenile salmon during emigration 
to the ocean. Currently, in-river survival during migration 
averages 40-60% per year for populations that emigrate as 
yearlings after one year of freshwater rearing, although certain 
conditions, such as very low river flow, can result in mean annual 
in-river survival around 25% [12,15]. However, overall survival to 
maturity remains relatively low (~ 1 %) for certain populations. 
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such as spring/ summer Chinook salmon, which return to their 
natal rivers in the spring and summer, from the Snake River, 
which is the largest tributary of the Columbia River. Other 

populations that also emigrate long distances (>500 river 
kilometers, rkm), such as the spring Chinook salmon from the 
Mid-Columbia River, experience greater survival (~3-4%) 
[12,13]. These consistent differences in survival led to a hypothesis 
that multiple stressors during migration in the Snake River 
increase juvenile mortality during estuarine and early ocean 
residence, a phenomenon referred to as "delayed mortality" 
[14,15]. However, it is not yet clear if yearling emigrants from the 
Snake and Mid-Upper Columbia River, which emigrate at 
comparable times and sizes, respond similarly during early marine 
residence, which is considered a critical period for anadromous 
salmonids [16,17]. 

Columbia and Snake River Chinook salmon that emigrate as 
yearlings enter the northern Cahfornia Current, an eastern 
boundary current characterized by seasonal upwelling. For the 
Mid-Upper Columljia River spring Chinook salmon genetic stock 
group, there- is c'\idence in support of the "growth-mortality 
hypothesis", which postulates that larger, faster growing individ- 
uals survive better due to their enhanced ability to capture prey 
and avoid predation compared with smaller counterparts [18,19], 
during their early marine residence. Tomaro et al. [20] 
determined that interannual variation in early marine growth 
and juvenile size after ~30 d in the coastal ocean was highly 
correlated with subsequent returns of Mid-Upper Columbia River 
spring Chinook salmon. However, there was weaker support for 
the "match-mismatch hypothesis", which posits that survival 
during early life is related to an appropriate overlap between 
predators and their prey [21,22]. 

For spring/summer (sp/su) Snake River Chinook salmon, 
which also emigrate as yearlings, prior research indicates that in- 
river juvenile survival is positively related to body size [23,24] and 
that earlier emigrants tend to survive at higher rates than later 
emigrants [25]. However, information on mechanisms of mortality 
during early marine residence is more limited and there is a lack of 
information on the relative importance of freshwater, estuarine, 
and marine factors on the survival of Snake River sp/ su Chinook 
salmon. 

We combined ocean collections of juvenile salmon with genetic 
stock identification and otolith structural and chemical analyses to 
evaluate the growth-mortality and the match-mismatch hypothe- 
ses during initial marine residence in Snake River sp/ su Chinook 
salmon. We also evaluated the ability of juvenile attributes and 
environmental indice^s during (^arly marine residence to account 
for interannual variation in survival. First, we characterized 
interannual variation in juvenile salmon attributes, such as size at 
marine entry and marine growth rate, and compared these 
attributes with estimates of survival. Second, we evaluated the 
relative support for the growth-mortality and match-mismatch 
hypotheses during early marine residence by comparing environ- 
mental variables, juvenile attributes, and survival. Third, we 
determined which environmental variables accounted for the 
observed interannual variation in juvenile attributes, specifically 
mass at capture. Finally, we used multiple linear regression and 
multi-model inference to evaluate the relative importance of 
biological and physical indices in describing interannual variation 
in survival. We used a retrospective approach to hindcast survival 
for years in which we had ocean survey data (1998-2008) and a 
forecast approach where we used the top hindcast models to 
predict survival during recent emigration years (2009-201 1). 



Materials and Methods 

Study Organism 

Snake River sp/su Chinook salmon was listed as "threatened" 
under the ESA in 1992. The current management unit, identified 
as an Evolutionarily Significant Unit (ESU), includes all naturally 
spawned sp/ su Chinook salmon in the mainstem Snake River and 
the Tucannon, Grande Ronde, Imnaha, and Salmon river sub- 
basins as well as numerous artificial propagation units [26]. An 
estimated population of > 1,000,000 fish in the late 1800s declined 
to <5,000 in the 1990s [26]. Survival of Snake River sp/su 
Chinook salmon, as indicated by smolt-to-adult return ratios 
(SARs) which provide a measure of overall survival from the 
emigrating smolt stage to the returning adult stage, ranged from < 
0.005 to 0.032 (or <0.5 to 3.2%) from 1998 to 2008 (http://www. 
fpc.org/). 

We used SARs for Snake River sp/su Chinook salmon as the 
survival metric for comparison with juvenile attributes and 
environmental indices. Snake River sp/su Chinook salmon 
emigrate as yearlings and the majority (>80%) spend two years 
in the ocean prior to returning to spawn [26]. SAR estimates 
include all ages-at-maturity for a particular brood year. We used 
the SARs that were based on dctex tions of emigrating smolts and 
returning adults with Passive Integrated Transponder (PIT) tags at 
Lower Granite Dam (LCD) on the Snake River (http://www.fpc. 
org/). We used the composite estimate for wild Snake River sp/ su 
Chinook salmon without jacks, which are precocious males that 
return after one year in the ocean, in subsequent analyses [12] 
(Fig. 1). We selected this metric because: (1) all SAR estimates for 
this population w(;re highly correJated from 1998-2008 (r = 0.986 
for SARs with and without jack; r = 0.740 for SARs for wild and 
hatchery composite; and r = 0.978 for SARs from LGD to 
Bonneville Dam and from LGD-LGD; (2) the available estimates 
for the proportion of jacks within a brood year were not 
significantly correlated with SARs (r = 0.455, 2000-2010); and 
(3) we were interested in making inference about the survival of the 
naturally-spawned (presumably wild) portion of the population, 
which is of primary conservation concern. 

Juvenile collection and genetic stock identification 

Surveys occurred off the coast of Washington and Oregon 
during late May and late June from 1998 to 2008 (Fig. 2). A 
Nordic 264 rope trawl was towed in surface waters, and catches 
were standardized to density (catch-per-unit-effort, CPUE) of 
yearling Chinook salmon based on trawl width and the distance 
towed (fish km~ ). On board, fish were placed immediately on ice, 
identified, measured (fork length (FL), mm), and sacrificed (i.e., 
frozen). AU animal work was conducted according to relevant 
national guidelines. Fish were collected under ESA Section 10 
permit #1410-7A, which is the federal process for research 
conducted by the National Oceanic and Atmospheric Adminis- 
tration (NOAA) that involves ESA-listed species. During this study, 
NOAA collections of fishes did not undergo a separate review by 
an Institutional Animal Care and Use Committee. 

In the laboratory, fish were re-measured and weighed (±0.1 g) 
and fin clips and otoliths were collected. Fin clips were used to 
genotype juveniles at 13 microsatellite DNA loci [27] and assign 
individuals to stock group using a standardized genetic database 
[28,29]. Stock assignments were made with the program ONCOR 
[30] and the likelihood model of Rannala and Mountain [31]. For 
Columbia River Chinook salmon, the stocks that are identifiable 
using the standardized set of microsatellite loci are mostiy 
concordant with the eight ESUs in the basin [27,29]. However, 
translocations of hatchery populations in the region have been 
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Figure. 1. Smolt-to-adult return ratios (SARs). Estimates of SARs for populations of spring and summer Chinook salmon from the (a) Snake 
River, (b) Mid-Columbia River, and (c) Upper Columbia River. SARs are presented as percentages without the inclusion of precocious males (jacks) and 
were acquired from the Fish Passage Center (httpV/www.fpc.org/). The run timing (sp = spring, su = summer) and the hatchery or river of origin are 
also included. 

doi:l 0.1 371/journal.pone.009981 4.g001 



very extensive [26] and in recent genetic studies, the Mid and 
Upper Columbia Ri\ er spring-run ESUs have been combined into 
a single stock [20,27,29]. From 1998 to 2008, 755 individuals from 
the coastal surveys were classified to the Snake River sp/su 
Chinook salmon genetic stock group. However, the years 1998, 
2001, and 2005 were not included in the analysis of jm enile 
attributes due to low CPUE and therefore small sample size (n< 
10). Therefore, 732 juveniles coUected in 1999-2000, 2002-2004, 
2006-2008 were included in subsequent analyses. 

Juvenile migratory attributes: size at and timing of 
marine entry and early marine growtli 

A subsample of the juveniles was selected for otolith structural 
and chemical analyses to determine size at, and timing of, marine 
entry as well as marine growth and migration rates. The temporal 
(May vs. June) and spatial (across transects) (Fig. 2) distribution of 
the subsample was similar to the overall catch (x^<27.5, p>0.05 
for all years except 2006-2008). In the 2007 and 2008 subsample, 
there was a slight bias towards juveniles collected along the 
Columbia River transect in May (27.5<x^<29.5). In 2006, there 
was an over-representation of juveniles from a northern transect 
(LaPush) in May and an under-representation in June (Queets 
River) (x^= 121.8). Overall, however, we consider the subsample 
representative of the ocean catch of juvenile Snake River sp/su 
Chinook salmon. 



Sagittal otoliths were removed, cleaned, and polished using wet- 
or-dry paper (240-2500 grit) and lapping film (1-30 (Xm) to expose 
the dorsal-ventral growth axis using standard procedures for 
elemental analysis [32]. Otolith Sr and Ca were measured along 
the dorsal-ventral growth axis using a VG PQ, ExCell inductively 
coupled plasma mass spectrometer with a New Wave DUV193 
excimer laser. The laser was set at a pulse rate of 7 Hz and 
travelled across the sample at 5 |a,m s~' with a spot size of 30 or 
50 |J,m. Normalized ion ratios were converted to molar ratios using 
standard procedures [33,34]. Instrument precision (mean percent 
relative standard deviation) was 4.5"/o for Ca and 4.7"/o for Sr 
across all samples and days (n = 50) and accuracy for Sr:Ca was 
4% (n = 5) based on USGS MACS-1. 

Image analysis was combined with Sr:Ca data to d(^tcrmine 
otolith width (OW) at marine entry and to estimate the date of 
marine entry [35]. For each individual, the OW at the time of 
marine entry was determined by the initial and abrupt increase in 
otolith Sr:Ca, which indicates exit from freshwaters, prior to 
stabilizing at marine values [36,37]. We enumerated the 
increments deposited after the initial and abrupt increase in 
otolith Sr:Ca to determine residence in brackish/ocean (hereafter 
referred to as "marine") waters. To determine date of marine 
entry, the duration of marine residence was subtracted from the 
date of capture. Marine migration distance was conservatively 
estimated as the linear distance between the mouth of the 
Columbia River (N 46.253°, W 124.059°) and the capture station 



PLCS ONE I www.plosone.org 



3 



June 2014 I Volume 9 | Issue 6 | e99814 



Snake River Spring/Summer Chinook Salmon Survival 




125* W 124' W 123" W 



Figure 2. Map of study location, (a) Columbia River watershed with locations of the mainstem dams and gaging station referred to in text. BON = 
Bonneville Dam; TDA = The Dalles Dam; JDJ = John Day Dam; MCN = McNary Dam; ICH = Ice Harbor Dam; LMJ = Lower Monumental Dam; LGS = 
Little Goose Dam; LGR = Lower Granite Dam. (b) Transect and station locations for ocean collections used in this study. 
doi:1 0.1 371/journal.pone.009981 4.g002 



plus 32.1 km to account for travel through the estuary. We divided 
the migration distance (km) by the marine residence time (d) to 

calculate the mean migration rate (km d~') for each fish, which 
was converted to body lengths per second (bl s~ ) based on 
estimated size at marine entry. 

We developed a direct back-calculation model based on the fork 
length to otolith width relationship of yearling sp/su Chinook 
salmon from the interior Columbia Ri\ er basin that were collected 
from 1999-2008 (r^ = 0.82, n = 362, p<0.001). We determined 
size at marine entry using Eq. (1). 

In(FLM) = 1.126( + 0.028SE)-In(OWM) -3.690( ±0.210SE) (1) 

where FLm = fork length (mm) at marine entry, and OWm — 
otolith width (|j,m) at marine entry. Marine growth rates (% d""^, 
mm) were then determined by subtracting estimated size at marine 
entry from size at capture, dividing by size at marine entry, and 
multiplying by 100. 

River, estuaty, and ocean indices 

We compiled indicators of river, estuary, and ocean conditions 
during juvenile emigration for comparison with juvenile attributes and 
survival. Data on daily discharge in the lower river were obtained from 
the United States Geological Survey (Site 14246900 at 46°N, 123°W). 
We characterized attributes of the Columbia River plume, defined 



using a cutoff salinity of 28, with the output of simulation databases 
(www.stccmop.org/datamart/virtualcolumbiariver) [38], including 
plume size (area of the plume surface and volume of the 3D plume) 
and location (expressed in terms of coordinates of the centroid of the 
surface plume) [39,40]. We hypothesized that conditions during 
emigration would be the most relevant to survival variation but 
physical indices were averaged across seasons (January to March, April 
to June, etc.) to identify the most appropriate period. 

We examined two basin-scale environmental indices, the Pacific 
Decadal Oscillation (PDO) and the North Pacific Gyre Oscillation 
(NPGO), which are statistically independent modes of variation in 
ocean sea surface temperature (SST) and sea level height (SLH), 
respectively. The PDO is defined as the leading principal 
component of North Pacific monthly SST variability poleward 
of 20°N [41]. Negative values of the PDO indicate cooler SST and 
relatively high salmon production off the west coast of North 
America [42,43]. The NPGO is the second leading principal 
component of variability in North Pacific SLH and is correlated 
with salinity, nutrients, and [:hlorophyll values [44]. Monthly 
mean values for these indices were downloaded from http://jisao. 
washington.edu/pdo/PDO.latest and www.o3d.org/npgo/data/ 
NPGO.txt. Physical indices were averaged across seasons (January 
to March, April to June, etc.) to identify the most appropriate 
periods. 
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To provide an indication of interannual variation in ocean 
productivity, we also examined the Copepod Community Index 
(CCI). The CCI is a numerical representation of all copepod 
species that are present in more than 5% of the samples collected 
biweekly 9 km offshore of Newport, Oregon using a 50-cm 
diameter, 202-|J,m mesh ring net towed vertically from 5 m above 
the sea floor to the surface. The values are rotated Axis 1 scores of 
a non-metric multidimensional scaling ordination of species 
abundance by sample date from 1996 to 2010 [45,46]. During 
spring and summer, negative CCI values indicate the presence of a 
"northern communit)'", i.e., boreal, neritic taxa that are large and 
lipid-rich, whereas positive values indicate the presence of an 
"offshore or southern community" comprised of smaller, relatively 
Upid-poor species [47]. The CCI may be indicative of the 
nutritional quality of the food web supporting juvenile salmon and 
their prey. 

In the northern California Current, there is a marked spring 
transition with the initiation of seasonal upwelling [48]. The 
physical spring transition is defined as the first day that the 10-d 
average for upweUing indices is positive and sea level height is 
negative [49] . The biological spring transition is defined as the first 
day that the copepod community present offshore of Newport, 
Oregon is dominated by lipid-rich, boreal species [45,50]. 

Evaluation of the Match-Mismatch and Growth-Mortality 

Hypotheses 

The majority of juvenile salmon along the west coast of the US 
enter marine waters during the productive spring/ summer 
upweUing season [51]. Hence, both the timing and magnitude of 
productivity during the upweUing season could be important to the 
growth and survival of juvenile salmon. There is evidence that, 
within years, juveniles survive better when they emigrate earlier in 
the spring [25]. However, it is not clear if interannual variation in 
survival is related to timing of in-river transit or marine entry, and 
there may be an interaction between the timing of marine entry 
and the initiation of the upweUing season. Therefore, in addition 
to our examination of river, estuary, and ocean conditions during 
juvenile emigration, we also estimated the date of marine entry in 
relation to the physical and biological transitions in coastal waters. 
For each individual, we determined the number of days between 
their marine entry and the physical (MEpx) and biological spring 
(MEbt) transition dates for that emigration year (marine entry - 
MEpT or MEbt) and generated annual mean estimates for 
comparison with survival. The expectation is that if the timing of 
marine entry of those juveniles that survived initial marine 
residence was important for subsequent survival, there would be 
significant, positive relationships between SARs and MEp^ O 
MEbt- 

In Pacific salmon and other fishes, survival is often positively 
related to body size at some point in the early life history 
[18,19,24]. However, it is not clear when the size-survival 
advantage occurs, e.g., during in-river migration or early marine 
residence. By determining size at marine entry and early marine 
growth rates, we can examine the relative importance of 
interannual variation in entry size and early marine growth on 
survival. Therefore, we related SARs to juvenile size at marine 
entry, early marine growth, and size at capture using correlation 
analyses (1999-2000, 2002-2004, 2006-2008). 

For aU correlations, we adjusted degrees of freedom when 
determining significance for Pearson product-moment correlation 
coefficients (^P) if needed to account for auto-correlation within 
time series, as recommended by Pyper and Peterman [52]. Given 
the high rates of marking Snake River Chinook salmon hatchery 
fish by clipping adipose fins and injecting coded wire tags across 



the study period (>90%; Regional Mark Information System, 
www.rmpc.org), we also compared attributes of marked and 
unmarked fish to determine if there were significant differences 
between these two groups. 

Hindcasting survival 

Our intention was to evaluate the relative importance of local 
and regional indices in describing interannual variation in survival 
as weU as to determine if juvenile attributes after a hypothesized 
critical period provide additional or unique information on 
survival variation (1999-2000, 2002-2004, 2006-2008). First, for 
the years with adequate numbers of juvenUes, we compared SARs 
with size at, and timing of, marine entry, marine growth rate, and 
size at capture to evaluate their relative importance in relation to 
survival using correlation analyses. Second, we used a multiple 
linear regression aj)proach to determine which environmental 
indices accounted for variation in those juvenUe attributes that 
were found to be related to survival (SARs) in our correlation 
analyses. Third, we evaluated the relative ability of the environ- 
mental indices to hindcast survival directly, including all years 
(1998-2008). Variables for model inclusion were selected based on 
visual inspection and correlation analysis and transformed when 
necessary to meet parametric assumptions. Variables with cross- 
correlations >0.50 were not included in the same model. We 
calculated Akaike Information Criteria adjusted for smaU sample 
size (AICc) to evaluate models [53]. Normalized likeUhood values 
{Wi) were calculated using the foUowing: 

exp( — 0.5 * A,) 



^jexp(-0.5* A,) 

where w,- are Akaike weights for model i, the numerator is the 
likelihood for model and the denominator is the sum of the 
relative likelihoods for models 1 to r. A,represents the difference 
between the AIC,. of the best model and the others. We also 
calculated model-averaged parameter coefficients and variance 
using standardized coefficients to evaluate the relative importance 
of variables in the models using the foUowing: 



vim- 



var(ft|g,) + (ft-ftr 



where fij denotes the model-averaged coefiicient, w,- is the model 
weight, and JS,- is the coefficient for model /. The unconditional 
variance estimate var(;S) was determined using the weight (w,) and 
variance (var(^, | g,- )) of each model [53]. 

Results 

Juvenile migratory patterns 

Across all years, the annual mean emigration date occurred 
from 6 May to 18 May and all juveniles emigrated between 20 
April and 19 June with a shift toward later emigration in 2006- 
2008 (Fig. 3). The mean proportion of marked, and presumably 
hatchery, fish in ocean collections was 80.5% and ranged from 
67% in 2000 to 89% in 2004 and 2006 (Table 1). Fourteen of 
these marked fish had PIT tags [54], which provided data on 
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juvenile size at release and information on when and where tiiey 
were detected within the Columbia River hydropower system. 
These tagged fish provided an opportunity to compare our otolith- 
derived estimates for size at, and timing of, marine entry with 
available information (http://www.ptagis.org/) (Table SI in File 
SI). Two fish (14%) were estimated to be smaller at marine entry 
than at release (by ~6%) and these were transported by barge 
through the hydropower system quickly (£10 d). Additionally, all 
of our estimates of marine entry were within five days of the fish's 
last detection within the hydropower system. Overall, these data 
provide a qualitative indication of the ability to estimate size at, 
and timing of, emigration. 



A relatively high percentage (34%) of the juveniles displayed no 
evidence of marine residence, i.e., no elevated Sr:Ca in their 
otoliths, which indicates recent marine entry (<5 d) [55,56]. 
Therefore, individual residence in coastal waters prior to capture 
ranged from 1 to 81 d with an overall mean of 20 d (± 1 7.9 SD). 
For those individuals with evidence of marine residence in their 
otoliths, mean annual marine growth rate ranged from 0.47% d~' 
in 2002 to 0.83% d"' in 2000 (Table 1). Mean migration rate 
ranged from 0.1 1 to 1.77 bl s ' and tended to increase later in the 
year (Fig. 4) (mean = 0.40 bl s"'±0.21 SD and 0.64 bl s"'±0.32 
SD in May and June, respectively, t-test, P<0.05). 




Day of year of marine entry 



Figure 3. Timing of marine entry for juvenile Snake River sp/su Cliinook salmon. Percent frequency for the day of year of marine entry 
(bars) as estimated from otolith chemical and structural analyses (n = 230). Year and mean day of year of emigration are included on each graph. Black 
lines represent cumulative frequency. Dotted lines represent overall mean date of emigration (May 12=133). 
doi:1 0.1 371/journal.pone.009981 4.g003 
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Figure 4. Juvenile Snalte River sp/su Cliinool< salmon marine 
migration rates, (a) Estimated marine migration rates (bl-s^^) for 
juvenile Snake River sp/su Chinook salmon across all years, (b) 
Individual marine migration rates for all years (1999-2000, 2002-2004, 
2006-2008). Filled circles represent juveniles collected during May 
cruises (n = 70) and open circles represent juveniles collected during 
June cruises (n = 77). Shaded boxes indicate cruise dates. 
doi:10.1371/journal.pone.0099814.g004 

There were significant interactions between year and presump- 
tive origin (marked vs. unmarked) for most juvenile attributes, 
including size at marine entry rate (Fy i3o>2.7, PO.05), marine 
migration rate (Fy i3o>2.4, P<0.05), marine growth rate (Fy i3o> 
2.4, P<0.05) and size at capture (Fy yjs = 2.6, P<0.001) (Table S2 
in File SI). There was no interaction between origin and year for 
the date of marine entry, and unmarked fish entered an average of 
2 days later than marked fish (Fi y = 7.0, P<0.05). We observed no 
other consistent differences between origins. 

Match-Misnnatch and Growth-Mortality Hypotheses 

We observed minimal evidence in support of the match- 
mismatch hypothesis (Table 2). Juveniles consistently entered 
marine waters after the physical spring transition but there was no 
statistically significant relationship with the physical or biological 
transition (Fig. 5). However, there was a non-significant positive 
trend (r = 0.639) with higher SARs when fish emigrated later in the 
year relative to the biological transition (Fig. 5). 

There was stronger support for the growth-mortality hypothesis 
after marine entry (Table 2, Fig. 6). Mean size at marine entry 
displayed negative, non-significant trends with survival but length 
(Fig. 6b) and growth (Fig. 6c), measured an average of 20 d after 
marine entry, were strongly, positively related to survival (Fig. 6c, 
r>0.73). 

Given that fish mass at capture was the most informative 
juvenile attribute in relation to survival (Table 2), we evaluated the 
ability of physical and biological variables to hindcast interannual 
variation in juvenile size. This analysis could only be completed for 
years with adequate ocean collections of juveniles (1999-2000, 
2002-2004, 2006-2008). Columbia River plume area and volume 
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Table 2. Pearson's correlation coefficients for comparisons between biological and physical indices and smolt-to-adult return 
ratios (SAR). 





FL (mm) 


M (g) 


SAR 


1. May marine density (yearling km"^) 


0.759 


0.773 


0.770 


2. June marine density (yearling km"^) 


0.797 


0.833 


0.886 


3. iWarine growth rate (% d"', mm) 


0.778 


0.805 


0.727 


4. Size at marine entry 


-0.100 


-0.142 


-0.252 


5. Juvenile migration rate (bl sec"^) 


-0.315 


-0.253 


-0.070 


6. Columbia River floW4_7 (m^ s"^) 


0.311 


0.310 


0.514 


7. IMEbt 


0.510 


0.547 


0.639 


8. MEpT 


-0.157 


-0.129 


0.090 


9. PD07_9 


-0.853 


-0.874 


-0.918 


10. NPG04_5 


0.906 


0.893 


0.715 


1 1. Plume area4_7 (km^) 


0.894 


0.901 


0.935 


12. CCIs 


-0.868 


-0.860 


-0.799 


13. SAR 


0.883 


0.902 


1.000 



FL = fork length (mm). M = mass (g). MEbt is the annual mean day of marine entry in relation to the biological spring transition whereas MEpj is the annual mean day of 
marine entry in relation to the physical spring transition (see text for additional details). Subscripts indicate the months over which data were averaged (e.g., PD07_9 = 
mean PDO from July to September). Adjusted significant values ('^P<0.05) are indicated by bold letters, n = 8 for all comparisons. Variables were In-transformed (FL, M, 
and (4), (6), (7) and (12)) or square root transformed (SAR and (2)) to normalize distributions and homogenize variances. Years included are 1999, 2000, 2002-2004, 
2006-2008. 

doi:1 0.1 371 /journal.pone.009981 4.t002 



are highly correlated (r = 0.963); therefore, we used plume area in 
our models due to its slightly better relationship with juvenile mass 
at capture (r = 0.901 vs 0.871). The top five out of 15 possible 
models accounted for similar proportions of the variance in 
juvenile mass at capture (^0.74) (Table 3). In general, yearlings 
were heavier in years in which the plume was larger, the PDO 
index was more negative, the NPGO index was more positive, and 
the CCI was more negative, i.e., dominated by northern, boreal 
copepod species (Table 2). Given the family of models, the model 



that incorporated plume area during emigration (April through 
July) was 1.4 to 2.5 times more likely than models based on basin- 
scale indices (NPG04_5 and PD07_9) and >3.5 times more likely 
than the models with the CCI5. As all four of these variables were 
correlated (r^O.59) and thus not included in the same model, we 
compared the model-averaged coefficients, which indicated that 

PlArea4_7 was the most informative variable ;8/ + var = 0.359 ± 
0.086), followed by NPG04_6 (0.260±0.057), PD07_9 (- 
0.141±0.022), and the CCIb (-0.041 ±0.042). 



Physical Transition (r = -0.090) Biological Transition (r = 0.639) 

-40-20 0 20 40 60 80 100 120 -60-40-20 0 20 40 60 80 100 120 



High 2008 

2000 

1999 

2002 

2007 

2006 
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Survival 2003 



I 0 1 




^ — I 



I — n — I 



-Oh 
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□ 25%-75% 
I 1%-99% 



Figure 5. Timing of marine entry for Snai<e River sp/su Chinool< salmon in relation to ocean conditions. Mean date of marine entry in 
relation to the (a) physical and (b) biological spring transition. Years are ranked from the highest to lowest adult returns and labeled by juvenile 
emigration year. 

doi:1 0.1 371/journal.pone.009981 4.g005 
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Figure 6. Relationships between survival and juvenile salmon 
attributes. Relationship between smolt-to-adult return ratios (SAR) 
and mean annual size and growth characteristics for Snake River sp/su 
Chinook salmon. Back-transformed SAR versus mean (±SE) (a) juvenile 
size at marine entry, (b) size at capture, (c) and marine growth rate. 
doi:1 0.1 371/journal.pone.009981 4.g006 



Local and regional indices and survival 

Given the relative importance of plume area in accounting for 
variation in juvenile mass at capture, plume area was included rn 
our initial comparison of models to hindcast SARs, which limited 
the analysis to 1999-2008 because no plume simulations are 
available for 1998. However, plume size was not included in the 
top ten models, therefore model comparisons were completed 
without plume metrics across emigration years 1998 to 2008. 
Interannual variation in SARs was relatively well-described (R^> 
0.70) by physical and biological conditions during emigration. In 
general, SARs were higher when the PD07_9 and the CCIfj were 
more negative, the NPG04_6 was more positive, and CPUEg was 
greater. The two top hindcast models included the PD07_9, and 
they both were >6.7 times more likely given the data than any 
model that included CPUEg (Table 4). PD07_g was the most 

informative variable (^/ + var= -0.692±0.217), foUowed by CCIg 
(-0.115±0.039), CPUEg (0.059±0.007), and the NPG04_6 
(0.023±0.002). 

Snake River sp/ su Chinook salmon SARs have remained low in 
recent years. Therefore, we predicted SARs for the 2009-011 
emigration years (Table S3 in File SI) using the family of models 
based on physical and biological variables described above to 
evaluate the utility of the variables identified in our hindcast 
models. There was a growing divergence among model predica- 
tions across these recent years (Fig. 7). Models with only the 
PD07_9 displayed greater error in prediction of SARs (50^48% 
of observed) compared with the NPGO (68-212% of observed) or 
local biological indices (CPUE =86-166% and CCI =80-241% 
of observed). The NPG04_6, CPUEfj, and CCIg more successfully 
captured the recent variability in SARs (Fig. 7). Therefore, it 
appears that although the PDO index was remarkably informative 
in past years, other physical (NPGO) and biological indices (CCIg) 
and attributes of the juveniles (CPUE) provide a more accurate 
indication of survival variation in recent years. 

Discussion 

We observed strong evidence for the growth-mortahty hypoth- 
esis during early marine residence of Snake River sp/ su Chinook 
salmon, which is a fmding similar to the adjacent Mid-Upper 
Columbia River spring Chinook salmon genetic stock group [20]. 
In years during which adequate numbers of juveniles were 
collected for analysis. Snake River sp/ su Chinook salmon SARs 
were positively correlated with early marine growth and size at 
capture. However, interannual variation in size at marine 
(brackish/ ocean) entry was not significantly related to subsequent 



Table 3. Model comparisons for juvenile mass. 





Model 


RSS 


AIC^ 


A/ 






PIArea4_7 


0.050 


-28.672 


0.000 


0.398 


0.812 


NPG04_6 


0.054 


-28.049 


0.623 


0.291 


0.797 


PD07_9 


0.062 


-26.860 


1.812 


0.161 


0.765 


CCIe 


0.069 


-26.067 


2.605 


0.108 


0.740 


CCIe, {ca^f 


0.029 


-23.589 


5.083 


0.031 


0.890 



Results for models describing variation in juvenile fish mass (g) after initial marine residence for Snake River spring/summer Chinook salmon based on ocean conditions 
during juvenile emigration. CCIg^ Copepod Community Index in June; PIArea4_7= mean plume area from April to July; NPG04_6= mean value from April to June; 
PD07_9=: mean value from July to September, RSS- residual sum of squares, AIC^ =^ Akaike Information Criteria adjusted for small sample size. A, represents the 
difference between the AIC^ of the best model and the others, u', indicates the relative likelihood of the model given the data. Variables were transformed (logarithm or 
square root) to normalize distributions and homogenize variances. 
doi:1 0.1 371 /journal.pone.009981 4.t003 
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Table 4. Model results for salmon survival. 





Model 


RSS 


AIC^ 


A,- 






1. PD07_9, CCIe 


0.0024 


-78.148 


0.000 


0.430 


0.830 


2. PD07_9 


0.0039 


-78.054 


0.095 


0.410 


0.724 


3. CPUEs 


0.0055 


-74.235 


3.913 


0.061 


0.731 


4. PD07_9, (PD07_9)^ 


0.0035 


-73.802 


4.346 


0.049 


0.748 


5. NPG04_6, CPUEe 


0.0041 


-72.163 


5.985 


0.022 


0.707 


6. NPG04_6 


0.0067 


-71.953 


6.196 


0.019 


0.520 


7. ecu 


0.0088 


-69.014 


9.134 


0.004 


0.373 


8. ecu, (CCU)^ 


0.0055 


-68.872 


9.276 


0.004 


0.507 


9. CRFIoW4_7 


0.0121 


-65.495 


12.653 


0.001 


0.135 



Comparison of models describing variation in survival of Snake River sp/su Chinook salmon based on conditions during juvenile emigration {1998 to 2008). PD07_9 = 
mean value from July to September; CC\^= Copepod Community Index in June; CPU Eg = catch of yearling Chinook (fish km"^} in June; NPG04_6= mean value from 
April to June; and CRFIoW4_7= Columbia River flow from April to July. RSS = residual sum of squares, AIC^ = Akaike Information Criteria adjusted for small sample size. 
A, represents the difference between the AlCc of the best model and the others, ir, indicates the relative likelihood of the model given the data. Variables were 
transformed (logarithm or square root) to normalize distributions and homogenize variances. 
doi:1 0.1 371 /journal.pone.009981 4.t004 



survival which indicates that, for juveniles that survived their first 
month at sea, early growth was a more important determinant of 
their subsequent survival than their size at initial marine entry. 
Our data indicate that there are survival advantages associated 
with faster marine growth and larger body size attained during the 
first 3-4 weeks at sea, which are potentially related to reduced 
over-winter mortality [16]. It is important to note that our 
approach did not examine selective mortality during in-river 
migration. We did not evaluate whether juveniles that are larger 
when they initiate downstream migration survive better, although 
there is evidence that this is true [24,25]. Size-selective mortality 
in-river would result in more uniform sizes at marine entry 
followed by high variation in growth that could influence survival, 
potentially during the subsequent winter. Our data reveal such a 
pattern with relatively uniform sizes at marine entry (mean 
= 134 mm FL with coefficient of variation (CV) = 14%) and 
greater variation in early marine growth (CV = 46%). Size- 
selective mortality requires individual variation in size and there 
may be multiple periods in the life history when size-selective 
mortality can occur, such as during in-river migration and during 
the first ocean winter. 

There was relatively weak support for the match-mismatch 
hypothesis, which is also similar to our finding for the Mid-Upper 
Columbia River spring Chinook salmon genetic stock group [20]. 
However, 80% of our juveniles were of hatchery origin, and their 
timing of marine entry can be influenced by hatchery manage- 
ment practices, such as release timing and transport (barging) 
protocols, in addition to natural variation in migration behavior. 
In general, juveniles consistently emigrated after the physical 
spring transition but there was minimal evidence that survival was 
greater if marine entry occurred longer after either the biological 
or physical transition. The lack of a significant relationship 
between survival and marine entry across years does not mean 
that, across evolutionary time scales, emigration timing has no 
influence on survival. Rather, we interpret this result to indicate 
that, under current management practices, there is only weak 
evidence for a relationship between survival and annual mean time 
of marine entry for juveniles that survived their first month at sea. 
Furthermore, it is important to note that there could be intra- 
annual variation in survival related to emigration timing [25], and 
survival advantages associated with timing may have been more 



evident under more protracted juvenile emigration, such as 
occurred historically [57]. 

There were some intriguing differences between the analyses 
that included a reduced number of years due to inadequate 
numbers of ocean-caught juveniles (1999-2000, 2002-2004, 
2006-2008) and those analyses that encompassed all years (1998 
to 2008). For years in which juvenile attributes were examined, 
their marine growth, size at capture, and subsequent survival were 
well-described by a suite of variables. In these years (1999-2000, 
2002-2004, 2006-2008), variables indicative of conditions within 
the ocean basin (PDO and NPGO), local environment (plume area 
and copepod community composition), and juvenile abundance 
(CPUE) were all significantly correlated with survival (Table 2). 
However, when the years 1998, 2001, and 2005 were included in 
the analysis, the PDO was the most informative metric, and 
yielded a model that was >6 times more likely than other models, 
such as those with the NPGO or CPUE. These differences among 
models indicate that when juveniles can be examined after their 
initial 3-4 weeks in the ocean, certain attributes, such as size and 
growth, provide a strong indication of year-class strength. It is 
important to note that, witii the exception of 2005, the SARs 
associated with those low juvenile collection years were not the lowest 
in our time series (Table S3 in File SI). However, those low catch years 
were rather unique: 1998 and 2005 were considered the worst ocean 
conditions for salmon growth and survival (http://www.nwfsc.noaa. 
gov/research/divisions/fe/estuarine/oeip/g-forecast.cfm) and 2001 
had the lowest Columbia River flow and plume size. These 
observations indicate that the reduced statistical importance of the 
plume metrics when all years could be included in the analysis may be 
because certain plume conditions (i.e., large voliune/area) may be 
important but not wholly adequate for good juvenile salmon survival. 
Therefore, for years with moderate ocean conditions, plume conditions 
are related to survival. However, when there are very poor ocean 
conditions, such as 1998 and 2005, the relative importance of the 
plume in relation to survival is minimized. 

The ability to understand and predict fluctuations in marine fish 
populations has been a primary research focus for > 1 00 years 
[3,8]. The ultimate goal of identifying a suite of parameters that is 
relatively easy to measure and provides robust forecasts of 
abundance has proven elusive. Initial success in identifying likely 
mechanistic linkages, such as prey abundance and distribution 
[4,21,22], environmental thresholds or windows associated with 
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Figure 7. Model comparison for Snake River sp/su Chinook salmon survival. Observed and predicted smolt-to-adult return ratios (SAR) for 
Snake River sp/su Chinook salmon in emigration years 2009-2010 based on the top eight models presented in text. Mean prediction and the upper 
and lower 95% confidence intervals are presented for each model (see Table 4 for additional details). The "Actual" SAR values for 2009 and 2010 were 
obtained from http://www.fpc.org/. The 201 1 SARs was estimated based on the relationship between SARs and adult returns of spring/summer 
Chinook salmon to Lower Granite Dam at a -2-yr lag {r = 0.815, 1998-2010). 
doi:1 0.1 371/journal.pone.009981 4.g007 



high survival [5], or promising combinations of environmental 
correlates [49,58,59], often fail to result in viable predictive models 
for extended periods of time. However, the prospect of gaining a 
mechanistic understanding that forms the basis for more robust 
model development continues to motivate researchers. In this 
analysis, we did not seek to optimize our abihty to predict recent 
survival (2009-20 1 1 emigration years); rather we compared how 
well our top hindcast models performed in forecasting. The top 
hindcast model, based only on the PDO, substantially overesti- 
mated survived, potentially by >400% given the expected SARs 
for the 201 1 cohort, but the two models with variables indicative of 
lower trophic levels (CCI) and juvenile abundance (CPUE) yielded 
the most accurate predictions, yet stiU overestimated recent 
survival by >170%. The use of model-averaged predictions has 
received increased attention recently and may prove useful 
[53,60]. However, given the high weight of the PDO (-0.692) 
in our family of mod(;ls, a model-averaged estimate would still 
have substantially over-estimated SARs for emigration years 2010 
and likely for 2011. The fact that the best predictor of recent 
survival was the catch of juveniles in June (CPUE) demonstrates 
that the acquisition of biological information on a population after 
significant mortality events, or critical periods, may provide some 
of the most accurate indicators of changes in regulatory 
mechanisms. However, effectively integrating such information 
into management structures remains an important challenge. 

Burke et al. [61] examined the relationships between adtilt 
returns of sp/ su Chinook salmon to Ice Harbor Dam on the Snake 
River and a suite of physical and biological variables. They also 
found evidence for bottom-up, growth-mediated influences on 
survival and highlighted the importance of basin-scale indices, 
particularly the PDO, but they cautioned that such basin-scale 
relationships can be regime-dependent. The mechanisms regulat- 



ing productivity and abundance can vary across climate regimes 
and, thus, predictive population models may be "regime-specific" 
[62,63,64]. Interestingly, the recent low survival of the 2010 and 
201 1 Snake River sp/su Chinook salmon cohorts (2012 and 2013 
adult returns) occurred during a period of strongly negative PDO 
values. Therefore, even within regimes, the relationship between 
survival and basin-scale indicators can vary substantially and 
shotild be interpreted with caution. 

The interannual patterns in size at marine entry, early marine 
growth, migration rate, and size at capture that we observed for 
the Snake River sp/ su Chinook salmon yearlings are very similar 
to those observed for the Mid- Upper Columbia River spring 
Chinook salmon yearlings [20]. For example, annual mean size at 
marine entry ranged from 127—150 mm FL for the Snake River 
and from 126-156 mm FL for the Mid-Upper Columbia River 
stock. Mean size at capture was also similar (147-169 mm FL vs. 
147-179 mm FL, for the Snake- River and Mid-Upper Columbia 
River stock groups, respectively) [20]. Given that the only 
consistent difference was that the Snake River yearlings entered 
the ocean an average of 7-10 days later than Mid-Upper 
Columbia River yearlings, the small differences in mean size 
could be related to duration of marine residence. Other studies 
focused on these two populations also reported similarities: 
Rechisky et al. [13] compared early marine survival of yearling 
Snake River and yearling spring Chinook salmon from the Mid- 
Columbia River using acoustic tags in 2006, 2008, and 2009 and 
reported a high level of covariation in early marine survival 
between these two interior Columbia populations. They suggest 
that the cause of the consistently lower overall survival for Snake 
River sp/su Chinook salmon when compared to Mid-Columbia 
spring Chinook salmon may occur north of southern Vancouver 
Island, BC, which was the northern extent of their detection array. 
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The lower survival of the Snake River sp/su Chinook salmon 
compared with the Mid-Columbia River spring Chinook salmon 
may be, at least in part, due to conditions further north and later 
in the life history. However, it is notable the two models that most 
accurately predicted recent SARs for Snake River sp/ su Chinook 
salmon included indices of conditions within Northern California 
Current coastal waters, cohort abundance (CPUE) and the 
copepod community (CCI), which indicates that local conditions 
are important for understanding overall survival. Furthermore, the 
observation that juvenile CPUE in both May and June were 
positively correlated with SARs across years from 1999 to 2008 
(0.759 and 0.855, respectively) provides additional evidence that 
conditions at marine entry or during very early marine residence 
influence subsequent survival [8,9]. 

Organismal-level studies focused on changes in size, growth, or 
condition of individuals before and after critical periods in the life 
history can provide valuable insight into likely mechanisms of 
mortality [65,66,67]. Our approach, which combined field, 
genetic, and otolith-derived information, provided novel informa- 
tion on early marine residence in an ESA-protected population. 
We determined that juvenile abundance and size during early 
marine residence and local (plume area) and basin-scale (PDO) 
indicators were all good indicators of subsequent survival (r>0.85). 
In the absence of information on juvenile attributes, basin-scale 
indicators accounted for a lower but still substantial amount of the 
variation in survival (r>0.70). Although the low survi\ al of Snake 
River sp/su Chinook salmon population may be related to factors 
within the river system and/or events that occur later in the life 
history, indices of cohort abundance and the copepod community 
within coastal waters remained the most informative of the 
available indicators in recent years. Future efforts to gain a 
mechanistic understanding of the population productivity of 
anadromous fishes will continue to benefit from organismal-level 
explorations across the life history. 
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